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Recently it has been argued that a three-body colour confinement interaction 
£T) • can affect the stability condition of a three-quark system and the spectrum of 

a tetraquark described by any constituent quark model. Here we discuss the 
[ role of a three-body colour confinement interaction in a simple quark model 

CN ] and present some of its implications for the spectra of baryons, tetraquarks 

and six-quark systems. 
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I. INTRODUCTION 



o In the strong coupling limit the SU<7(3) colour gauge group leads to a three-body colour 

confinement interaction in baryons In SU(3) lattice QCD the static three-quark potential 
can be measured with the help of 3-quark Wilson loop operators (see for example 
In recent lattice calculations the ground state potential of a three-quark system has 

been extracted as a sum of a two-body Coulomb plus a three-body interaction confinement. 
These studies lead however to rather contradictory results; one of them gives support to 
the Y-type flux tube picture of Ref. [|IJ, while the others favour the A ansatz (where the 
three-quark potential consists in a sum of two-body components). 

The interaction potential obtained in these calculations corresponds in any case to the 
colourless ground state only and no information from lattice QCD about colour octets is 
available so far. In practice, for simplicity, in quark models as e.g. that of the confinement 
is treated approximately as a two-body colour operator. This can be expressed in terms of 
the quadratic (Casimir) invariant operator of SU(3). 

This three-body colour confinement interaction should not be confused with the three- 
body force |7]|| associated to the instanton 't Hooft's interaction, which in the nonrelativistic 
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limit contains a colour operator similar to the one introduced below, but is of short range, 
in contradistinction to confinement forces. The instanton-induced three-body force cancels 
for three quarks in a colour singlet state and is moreover only effective if the three particles 
are in a flavour singlet state. Another short-range three-body force has also been introduced 
in || on a purely phenomenological ground. This force has a simple scalar structure and 
accounts for a better description of the Roper resonance. 

Based on the algebraic argument that SUc(3) is an exact symmetry of QCD, which 
implies that any quark model Hamiltonian inspired by QCD can be written in terms of SU(3) 
invariant operators, a three-quark potential that depends on the cubic invariant operator of 
SU(3) has recently [0J been added to the usual two-body confinement. Its implications on 
the spectrum of ordinary g 3 and exotic q 2 q 2 hadrons have been considered. In particular 
turning on a three-body force with an appropriate strength denoted by c, it was shown that: 
1) the three g 3 colour states, namely 1, 8 and 10 appear in the correct order, 1 being the 
lowest one, as it should be, 2) in the q 2 q 2 system the three-body interaction brings distinct 
contributions to the two possible colour singlet states, by enhancing the binding in one and 
diminishing it in the other, depending on the sign of c. 

In this study we rederive some of the relations found in Ref. [[UJ and discuss explicitly 
the role of confining three-body forces in g 3 , q 2 q 2 and g 6 systems. In the simple framework 
of a harmonic confinement we show that in the g 3 system there is a competition between 
the three-body force and the kinetic energy in rising the energy of the colour octet states. 
By introducing singlet- singlet and octet-octet gg coupled pairs we show that the range of 
values of the strength of the 3-body force giving a correct spectrum for a g 3 system also 
favourably affects the spectrum of a tetraquark q 2 q 2 . The results on the role of a 3-body 
confining force in a g 6 system, relevant for the NN problem are entirely new. 



II. THE THREE-BODY FORCE FOR BARYONS 



In this section we recall and discuss the findings of Ref. |10[ in relation with a three-body 
interaction of type 

V36 = Vijk = VijkCijk , (1) 

with 

V ijk = \cmu 2 [(n - r 3 ) 2 + (r, - r k f + (r k - nf] , (2) 

where c is a strength parameter and Cij k a colour operator of type 

C ijk = d abc F? F* F c k , (3) 



i j fc ' 
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where F® = -A" is the colour charge operator of the quark i and d a c some real constants, 

symmetric under any permutation of indices and defined by the anticommutator of the 
Gell-Mann matrices A a as 



{A a , X b } = 2d abc \ c . (4) 
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These constants satisfy the following orthogonality relation 

d abc d abe = 5 ^ (5) 

The operator ([3]) can be expressed in terms of the two independent invariant operators of 
SU(3) as [0 (f° r a proof see Appendix A) 



d a bc F a F > ps= l [ C (3) ^ _ 5^ + », (6) 



where we slightly changed the notation of jnj, by writing C^ 2 -* instead of for the 
quadratic invariant and instead of for the cubic invariant. For a given irrep of 
SU(3) labelled by (A//), the eigenvalues of these invariants are 

(C&) = i(A 2 + fi 2 + Xfi + 3A + 3/i) (7) 



and (see for example Refs. fl ij or jTJ] ) 

(C (3) > = 1(A - a*)(2A + a* + 3)(A + 2// + 3) . (8) 

Then for a g 3 system the expectation values of (H) are — — , - for a singlet (A//) = (00), 

9 3o 9 

octet (A/i) = (11) and decuplet (A/i) = (30) states respectively These are the coefficients 
appearing in the last term of Eq. ( |14D below. 

Turning on a 2-body confining interaction, which ensures stability for a qq pair and 
adding the 3-body confining interaction (E)-(H) of strength c relative to the 2-body one, in 
Ref. [|10J it was found that the spectrum of a q 3 system is correctly described provided 

3 2 

--<c<- (9) 

The closer c is to the lower limit, the larger is the gap between the colour octet and singlet 
states. To see this, let us consider the Hamiltonian 

H = T + V 2b + V 3b (10) 

where T is the kinetic energy and V^b a 2-body confinement interaction of the form 

V 2b = J2V ij ( Cl + t.+F l a F°) (11) 
i<j 6 

containing an arbitrary constant c\ which we set equal to 1 as in Ref. and take 

Vij = - mu 2 (n - rj) 2 (12) 

V36 is the 3-body confinement interaction of Eqs. ([l|)-(0). Performing integration in the 
colour space and expressing H in terms of the internal coordinates p — (f{ — r2)/\/2 and 
A = (fi + f 2 — 2f* 3 ) I Vq, we have: 
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H = 3m- A-(VJ + V 2 ) + \mu 2 Xl (p 2 + A 2 ) 



with 



(13) 



X, 



= < 
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i=l (singlet) 
i=8 (octet) 
i=10 (decuplet) 



(14) 



In the expressions of \% (i — 1>8 or 10), the first and second terms stem from the colour 
part of V26 and V3& respectively. 

We search now for solutions of H. In order to satisfy the Pauli principle, in the lowest 
colour-singlet state the quarks should be in a s 3 configuration. This implies that the orbital 
wave function is symmetric and one can take: 



1 



'00 



^3/253 



exp [-(p 2 + A 2 )/(26 2 )] 



(15) 



with b a variational parameter. 
The expectation value of H is then: 



The minimization with respect to b 2 gives: 



(16) 



(17) 



with bl 



h 



mw 



. The energy of the singlet state is thus: 



E l = 3m + 



3^ 
mb\ 



3m + 3huJ3xi 



(18) 



The colour-octet 8 must be combined with an s 2 p configuration in order to satisfy the 
Pauli principle. Indeed, the spin-flavour part being totally symmetric, as for the nucleon 
ground state, the orbital-colour part must be antisymmetric, i.e. it is of the form: 



where 



0?o 



^0 



7T 3 6 8 

2 

7T 3 6 8 



1/2 



1/2 



p z exp [-(p 2 + A 2 )/(26 2 
A, exp [-(p 2 + A 2 )/(26 2 



(19) 

(20) 
(21) 
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are the mixed orbital symmetry states with one unit of angular excitation and C p , C x their 
colour counter-parts. The corresponding eigenvalue of H is: 

2h 2 

E 8 = 3m + — + Qmuj 2 b 2 X8 (22) 
mb z 



and the minimization with respect to b 2 gives: 

_b 



bl = 4= (23) 



which leads to 



Ah 2 



E s = 3m + — 2 = 3m + Ahu J3xs (24) 

7770 q 



Accordingly the gap AE between the octet and singlet states is: 



AE = Mty3 X 8 - 3 V 3 ^) ( 25 ) 

The largest value of AE corresponds to Xi — 0, i.e. to c = —1.5, in which case AE ~ 10.7hu. 
For c = (no three-body force) we would have AE ~ 3.5hu. This means that the gap is 
enlarged by a negative c and triples for the limiting value c = —1.5. 

The colour state 10 requires an antisymmetric orbital state in order to satisfy the Pauli 
principle. Its form is (see e.g. Ref. |12[], chap. 10): 

K = ) V2 (P + A- - P-X+) exp [-(p 2 + X 2 )/(2b 2 )} (26) 

with p± = p x ±ip y , etc. The subscript 10 means total L=l, M=0 as above. This is the only 
value of L allowed by an antisymmetric state built from the configuration sp 2 . Proceeding 
in a similar way as for the two previous cases, one gets: 

6 -7fc (27) 

and hence 



E w = 3m + 5huy/3xio (28) 
The gap between the decuplet and the octet state is thus: 



E 10 -E 8 = nu{hp X io - 4^3x8) (29) 

For c = —1.5, one has E w — E 8 — 3hu>, i.e. this state is located above the octet, as expected, 
with quite a large gap for a limiting value of c. 

Let us now evaluate the octet-singlet gap (p5l ) by taking c = -1.43 as in fK|. There is 
of course some arbitrariness in choosing hu. One can take for example b = 0.437 fm and 
m = 0.340 GeV, as typical values for a quark model (see e.g. Ref. |15|]), which give %u = 0.6 
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GeV. This implies a gap AE ~ 5.5 GeV. For c = (no three-body force) one would have 
AE = 3.5hu ss 2.1 GeV. 

Before ending this section we should note that there is some arbitrariness in fixing the 
lower limit of the range of c as given by the inequality @. This limit is related to the 
choice of the arbitrary constant c\ which has been set equal to 1 in Eq. (|XT|) . But taking 
for example c\ = 4/3, which is another good choice in constituent quark models, one gets 

Xi = 2 H c. The stability condition for the singlet would then gives 

9 

- I < c (30) 

i. e. a different lower limit, slightly more favourable than the one of inequality flj^) because 
with c = —9/5 and the new expressions for Xi, Xs one gets AE ~ 11. 5huj. 



III. THE THREE-BODY FORCE FOR TETRAQUARKS 

If F a is the colour charge operator of a quark, for an antiquark we must have 

F a = ~X a * (31) 

in order that F a (a = 1,2,. ..,8) satisfy the Lie algebra too. Then one can write the three-body 
interaction acting in the q 2 q subsystem as 

C ijk = -d abc Ft Fj T k (32) 

and the three-body interaction acting in the qq 2 subsystem as 

C ijk = d abc F? F) T k (33) 

As a remark, the three-body interaction in an antibaryon should be 

C ljk = -d abc F- F) T k . (34) 

In Ref. |H| the operator (|32|) is given in terms of SU(3) invariants as 



where acts on the q 2 q subsystem but the Casimir operator acts only on the subsystem 

of i+j quarks (see Appendix A). If the quark subsystem is in a symmetric state it gives rise 
to a q 2 q [211]c state called s and if it is in an antisymmetric state to a [211]c state called 

a. Both these states have (A//) = (10), which according to (H) gives (C^ r ) = — . But 

for the subsystem of the i + j quarks only, the SU(3) representations are different. One has 

(A/i) = (20) for the s state and (A/i) = (01) for the a state. Then the expectation value of 

5 5 

the operator ( |35"1) is — for s and — for a, consistent with Table II of HlO . 
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The operator fl33|) acting on a qq subsystem can be brought to a form similar to (p5|). 
This is 

The difference with respect to a q 2 q subsystem is now that one has to calculate the expec- 
tation value of C j+ j + ^ for a [221]^ colour state for which (A/i) = (01) so that one has now 

((7i+l+fc) = . The subsystem of antiquarks gives for C (2) the same value as that for the 



quarks so that the operators (^) and (33) have the same expectation value which leads to 
Eqs. (26) and (27) of Ref. [[TtJ 

Vs = Ys c mu}2 ^ 2 + r * 3 + r * 4 + r ™ + r24 + r ^ ^ 

V a = --c muj 2 (r 2 2 + r 2 3 + r 2 4 + r 2 3 + r 2 4 + r 2 4 ) (38) 

One can see that the contribution of the fourth particle is also included in these equation. 
To understand this one can for example add another antiquark to q 2 q. This leads to the 
singlet colour state [222] appearing from the direct product [211] x [11]. By construction 
this singlet has an intermediate coupling both between quarks and antiquarks. The two 
quarks couple either to a 3 or a 6 state and the antiquarks to 3 or a 6 state. If the particles 
1 and 2 are quarks and 3 and 4 are antiquarks, V a and V s are the three-body contribution 
to the colourless states denoted by |3i 2 3 34 > and (612634 > respectively. 

For a negative c, as required for baryons described by a pure constituent quark model (no 
gluon components in the wave function), the mass of the 1 612634 > state is reduced and the 
mass of |3i2334 > is enhanced by a three-body force. The situation is opposite for a positive 
c. In rpTof] , V s was associated to the unobserved sextet-sextet state, which would mean that a 
positive c is preferable. The conflict can be solved by noting that the relevant states in the 
present problem are in fact linear combinations of |3i 2 3 3 4 > and 1 612634 >. Such states are 
important asymptotically and they are defined by the transformations (see e.g. Ref. [|i"4| ) 




I13I24 >= \ o I 3 i23 3 4 > +\ ~ I612634 > , (39) 




13824 >= -\/^|3i 2 3 3 4 > +W^|6i 2 6 34 > ■ (40) 



In these states the intermediate coupling in the colour space takes place between a quark q 
and an antiquark q. This gives colour singlet qq pairs in Eq. ( [39|) and colour octet ones in 
Eq. fl40D . Asymptotically the energy of 1 813824 > must become large, as such a state is not 
expected to be seen. Using the transformations (|39| ) and ( [40]) one obtains the contribution 
of the three-body interaction in a q 2 q 2 system as 



(I13I24IC123II13I24) cc + ~ A] c = (41) 
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and 

(8 13 8 24 |Cl 2 3|8l3824> oc [-(-5) + I 1] c = - A c (42) 

which shows that with a negative c one raises the expectation value of the octet-octet above 
the singlet-singlet state, more than with c = 0. This implies that the coupling between 
octet-octet and singlet- singlet states due to a hyperfme splitting will be diminished, which 
amounts to make a ground state tetraquark less stable. This seems to be consistent with 
the experimental observation that no stable tetraquark system has been seen so far. 



IV. THE NN INTERACTION 

The short-range NN interaction can be studied as a q 6 problem. First we give a simplified 
discussion by considering that the six quarks are in a totally symmetric orbital state [6]o- In 
such a case the spin-flavour part of the wave function has a [33] fs symmetry which combined 
with the colour symmetry [222] q state leads to a totally antisymmetric state. The latter is 
a superposition of five colour components given by the five following Young tableaux: 





1 


4 




1 


3 




1 


3 




1 


2 




1 


2 




2 


5 


, ^2 = 


2 


5 


, ^3 = 


2 


4 




3 


5 


, ^5 = 


3 


4 




3 


6 




4 


6 




5 


6 




4 


6 




5 


6 



(43) 



Below we give some details of our calculations of the three-body matrix elements for a 6q 
system. In a state of orbital symmetry [6]o, i-e. of configuration s 6 , all orbital matrix 
elements are equal so one has only to calculate the colour matrix element: 

\ E E E < Ud ahc F^F^ n > (44) 

i<j<k a,b,c n=l 



where the factor 1/5 comes from the normalization of the total wave function (see e.g. [12 
, chapter 10). In Appendix B we have explicitly proved that: 

J2 < fa\d abc F?F*F°\il) 2 >= -5/36 (45) 

a,b,c 

We get the same result for ip 3 , ?/> 4 and ip 5 but not for ipi, for which one has: 

< ihld^F^F^lipx >= 10/9 (46) 

a,b,c 

i.e., the result for the singlet (123), as expected. As shown in Table 1, the matrix elements 
< ^ n \d ahc F^F^F^ n > differ for different (ijk) but the sum over the five states doesn't 
depend on the choice of (ijk); the table has been calculated for a given value of the colour- 
indices (abc) = (146) but the conclusion is also true for the other values of the colour 
indices. Then the calculation of the matrix element ([44]) reduces to the calculation of the 
matrix element of the three-quarks (123) which has to be multiplied by Cq = 20. 
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Thus we obtain 

i 6 5 2D 5 

5 E E E < ^n\d abc F^F^ n >= T E E < i>n\d abc F?F*F^ n > 

i<j<k a,b,cn=l a,b,cn=l 

= 4{J2< Hd^FfF^F^ > +4 < Hd abc F?F%F^ 2 >) ^ 

a,b,c a,b,c 

= 4(10/9 - 20/36) = 20/9 

Therefore we can see that in the case of a three-body confining force the value of the matrix 
element of a six-quark system in the symmetry state [6]o is equal to two times the value 
for a single baryon. This situation is similar to the two-body force where the expectation 
value of a 2-body operator V 2b cx 2J F?F? in the symmetry state [6]o is equal to - 4 i.e. two 

i<j 

times the value of a single baryon. Let us define an adiabatic NN potential as the difference 
between the interaction Hamiltonian at zero separation distance and at infinity, i.e. 

V NN = H(0) - H(oo) (48) 

In this difference only the kinetic energy survives if the NN system is in the state [6}o- The 
contribution of the confinement due both to two- and three-body forces cancels out because 

V conf (0) = V conf (oo) = V 2b + V 3b = -4 + 20/9c (49) 



The two-body confinement force has been discussed for example in Ref. |1| where the 
Hamiltonian also contains a hyperfine interaction. 

However the physical NN state is a combination of three symmetry states given by [16| 

\NN) = ^ \[6}o[33) FS > +^| |[42] [33] FS > |[42] [51] F5 > (50) 

For SI=(01) or (10), one should also consider the physical AA state: 

/~4~ /To 7 l~25 

|AA) = J- \[Q}o[S3} FS > +J- \[42} [33} FS > +J- \[42]o[51] F8 > (51) 



The unphysical colour octet-octet (CC) state has the form JT6 



\CC) = ^ \[Q]o[S3] FS > -yj |[42] [33] F5 > (52) 

By using 3-body fractional parentage coefficients (cfp) (given in Appendix C), we 
calculated the expectation value of the 3-body potential acting on the symmetry states 
| [42]o[33]fs > and |[42]o[51]f5 >• For |[6]o[33]fs >, the result is straightforward as shown 
above. In short, we have found the following expectation values: 

20 

< [Q]o[33] FS \V 3b \[6]o[33] FS >= —c (53) 

< [A2}o[33\ FS \V 3b \[A2}o[33} FS >= (54) 

< [42]o[51] FS |V^|[42]o[51] FS >= (55) 



Note that the | [42]o[33]fs > and | [42] o [51] fs > states have the same expectation value, 
consistent with the fact that is spin-isospin independent. Using the transformations 
( |50| - |5"2| ) from the symmetry states to the physical states NN, AA and the hidden colour 
CC state, we obtain the following matrix for V 3 t : 





NN 


AA 


CC 


NN 


28 
— c 
81 


38\/5 

c 

405 


38\/5 

c 

135 


AA 


38^ 
405 ° 


121 
405° 


76 

135 C 


CC 


38^ 
135 


76 

135° 


9 

-c 

5 



(56) 



The eigenvalues of this matrix are E\ = c/9,E 2 = c/9 and E 3 = 20c/9. This shows that 
the effect of the 3-body colour confinement on NN and AA is identical and rather small as 
compared to that on CC. In particular for a negative value of c, the spectrum of NN, A A 
and CC lowers and shrinks. For a positive c, the situation is the other way round. This 
means that, for c < 0, itself brings some attraction and implies a stronger coupling of 
CC to NN and AA due to a hyperfine interaction. This will lead to a reduced hard core 
repulsion in the NN potential. 

V. CONCLUSIONS 

We discussed the role of a schematic three-body confinement force in the spectra of 
3q, q 2 q 2 and q 6 systems. We found that a three-body confinement interaction with a negative 
strength c has the following effects: 1) it increases the gap between the physical colour 
singlet state and the unphysical coloured states in baryons; 2) it raises the expectation value 
of qq pairs in colour octet-octet states with respect to singlet-singlet states in tetraquark 
systems, which will lead to a smaller binding in a ground state tetraquark when a hyperfine 
interaction is included and 3) it increases the coupling between physical states and CC states 
in q 6 systems. While, in the first two cases, the gap between the physical colour-singlet state 
and the non-physical coloured states is increased, the opposite is true for the 6q system. The 
larger coupling between physical and CC states induced by the 3-body interaction has both 
negative and positive consequences: it will reinforce the undesirable Van der Waals forces 
but, on the other hand, it brings more attraction into the NN potential. This may be a 
desired feature for quark models which give a too strong hard core repulsion. 

A three-body force with a negative strength c will have just opposite effects than the 
ones mentioned above. 

Our results are valid for any quark model, irrespective of the hyperfine interaction. It 
would be useful to extend this study to a more realistic confinement interaction. 

Note also that our conclusions for q 2 q 2 and q 6 correspond to a zero separation between 
the interacting clusters (2 mesons and 2 baryons respectively). It may be possible that the 
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contribution of the confinement interaction changes with the separation distance. This is 
the aim of a further study. 



APPENDIX A: 



In this Appendix we first prove Eq. (|6]). In order to avoid any confusion, we fix the 
indices (ijk) of @ to be (123). We rewrite Eq. (|3]) as 



1 3 

d abc F a F b F c = -[J2 d abc F?F*F£ 



-3 dabc F i F i F j + 2 E dahC F i F i F i\ ( Al ) 

i,j i 

where the second sum in the right-hand side compensates for the extra terms contained in 
the first sum, but as we extract too many we add the third term for satisfying the equality 
correctly. The first term is precisely the cubic invariant operator acting on the three-quark 
system C\^_ 2 +z an d the last term is 2 x 3 times the cubic invariant operator C[ acting on 

a quark. The latter is replaced by its eigenvalue — so we get 

dabc papbpc = 1 [C (S) _ 3 £ d abc papbpc + 20 (A2) 
6 i,i 3 

Due to the fact that the constants d abc are symmetric under the permutation of indices we 
can modify the second term as 

jabc rpa rpb rpc \ 

« ^1-^2-^3 - g [°l+2+3 

3 3 90 
-gE^W.^+yl (A3) 

and simplify it by using the anticommutator (§) in the form 

W *?} = ^ 6c ^ c (A4) 
and the orthogonality relation (||). This leads to 

d abc F i F 2 F 3 = g [ Ci+2+3 _ 2^i+ 2 + 3 + "3" ] ( A5 ) 

i.e. Eq. (|). 

Next we prove Eq. (|35|) . One can rewrite the operator ([32|) as 
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- d abc F«F b 2 Fl = -i [d abc (F? + F* + F a 3 )(F b + F b 2 + f\){F^ + F 2 C + F%) 

2 2 

-3 J2 dabc F i F i F j + 2 E dabc K ap i F i 

-sj:d abc F?F b F c 3 

i 

-sj2d abc F?F b 3 F c 3 

i 

-d°» c Fp*F c 3 ] (A6) 

The operator in the first term of the right-hand side is the C^> invariant associated with 
the whole system formed of the quarks 1 and 2 and the antiquark 3. As in Eq. ( |A1| ), 
the extra-terms introduced by this operator must be compensated in order to recover the 
left-hand side. This is the role of the other terms. 

The first term can be replaced by C { +2 + 3 anc ^ ^ e secon d by 5/2 C} +2 where the factor 
5/2 has the same explanation as in Eq. (A5|). The third term is 2 x 2 the invariant operator 
for a single quark = 10/9. The last term contains only antiquark charge operators and 
is thus identical to = - 10/9. The sum over the constant terms gives 

4C® - Cf = 50/9 (A7) 
where q = 1 or 2 and q = 3, so we have 



d a bc F a F , F c = _ 1 [c (3) 2+ _ _ | c (2) 2 _ 3 £ ^ _ 3 £ ^ + ^ ( ^ 



Due to the fact that cf bc are symmetric under permutation of indices a, 6 and c we can again 
use the identity 

d abc F?F b F c 3 = \ d abc {Ft, F b }T 3 (A9) 

and 

d ahc F*F b 3 F c 3 = \ d abc F?{F*, F° 3 } (A10) 
From (|) and (^TJ) it follows that 

(F a l ,F b } = -d abc F c l . (All) 

From (K4) and (|A11| ) it follows that the third and fourth term compensate each other and 
we get 

^s = -i[C« +3 -^a + f] (A12) 

i.e. equation (|35|). 
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APPENDIX B: 



In this Appendix we give the details leading to Eq. (BSD. We first present the explicit 
expressions for the colour wave functions corresponding to the Young tableaux of Eq. (pL3"D . 
In the state ipi, the sets of particles (123) and (456) are both in a totally antisymmetric 
state. Therefore, one can write: 



r(l) 6(1) g(l) 
r(2) 6(2) g(2) 
r(3) 6(3) 5 (3) 



r(4) 6(4) 5 (4) 
r(5) 6(5) g(5) 
r(6) 6(6) 5 (6) 



(Bl) 



where r, 6, g denotes the different quark colours. By applying the permutation (34) to ipi, 



one gets 12 



from where one obtains : 



(34)^1 = + ^pV>2 



(B2) 



'02 



2V2 



4V2 



(34)0! - -Vi 



r(l) 6(1) 5 (1) 
r(2) 6(2) 5 (2) 
r(4) 6(4) 5(4) 



r(3) 6(3) 5(3) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(3) 6(3) 5(3) 



r(4) 6(4) 5(4) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



(B3) 
(B4) 



The states 1^3,^4 and ^5 can be obtained by a similar procedure. They read: 



2VQ 



r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(5) 6(5) 5(5) 



r(3) 6(3) 5(3) 
r(4) 6(4) 5(4) 
r(6) 6(6) 5(6) 



r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(4) 6(4) 5(4) 



+ - 



r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(3) 6(3) 5(3) 



r(4) 6(4) 5(4) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



r(3) 6(3) 5(3) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



(B5) 



2v/6 



r(l) 6(1) 5(1) 
r(3) 6(3) 5(3) 
r(4) 6(4) 5(4) 



r(2) 6(2) 5(2) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(4) 6(4) 5(4) 



1 

+ 2 



r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(3) 6(3) 5(3) 



r(4) 6(4) 5(4) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



r(3) 6(3) 5(3) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



(B6) 
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^5 



3V2 



r(l) 6(1) 5 (1) 
r(3) 6(3) 5 (3) 
r(5) 6(5) 5 (5) 



r(2) 6(2) 5 (2) 
r(4) 6(4) 5(4) 
r(6) 6(6) 5(6) 



r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(5) 6(5) 5(5) 



r(3) 6(3) 5(3) 
r(4) 6(4) 5(4) 
r(6) 6(6) 5(6) 



r(l) 6(1) 5(1) 
r(3) 6(3) 5(3) 
r(4) 6(4) 5(4) 



r(2) 6(2) 5(2) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 

r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(3) 6(3) 5(3) 



1 

+ 4 



r(l) 6(1) 5(1) 
r(2) 6(2) 5(2) 
r(4) 6(4) 5(4) 



r(3) 6(3) 5(3) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



r(4) 6(4) 5(4) 
r(5) 6(5) 5(5) 
r(6) 6(6) 5(6) 



(B7) 



We want now to calculate the expression 

J] < Hd abc F?F$F£\ip 2 >. This can be done 



a,b,c 



explicitly using Eq. ( |B4|) . The non vanishing terms in this sum are given in Table [TI|. Taking 
into account the multiplicity of each term, one can directly check that the final result is -5/36. 
Similar calculations can be done for the functions ^3, ^4 and ^5 leading to the same answer. 



APPENDIX C: 



In this appendix, we give the values of the 3-body coefficients of fractional parentage 
(cfp) necessary to calculate the expectation value of the 3-body potential and we sketch the 
method to determine them. 

The state | [42] [33]^s > can be decomposed as: 



|[42] G [33] F5 >= 



1 


2 


.3 


4 


5 


6 



1 


3 


5 


2 


4 


6 



oc 



FS 




1 


2 


5 


3 


4 


6 



oc 



FS 




1 


3 


2 


5 


4 


6 



1 


2 


4 


3 


5 


6 



oc 



FS 




1 


2 


4 


3 


5 


6 



oc 



FS 



1 


4 


2 


5 


3 


6 



1 


2 


3 


4 


5 


6 



oc 



FS 



(CI) 



One has to determine the 3-body cfp associated to the decomposition of the OC part of the 
wave function into its orbital and colour parts, for example: 









4 


5 


6 







5 


6 




4 






oc 









4 


5 


6 



(C2) 



c 



To determine the 3-body cfp we need to write the Clebsch-Gordan (CG) coefficients of 
Sq specifying the place of the last three particles (pqr), where p,q,r represent the row in the 
Young tableau where the particles 6, 5 and 4 are located. The position of the remaining 
particles is denoted shortly by y. By using the factorization properties of the CG |L7[], one 
gets the following relations : 
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S ( Wq'r'y' [f"] P "q"r"y" \ [f]pqry) = K([f]p> [f]p" \ [f]p) 

x5([/;,]g'rV[^,]g'VV|[/ p ]gr 2/ ) 

x5([/; v ]rV[^ g »]rV|[/ P #Z/) (C3) 
= K([f]p'[f"}p"\[f]p) K {[QqV^q'VM) 

X S (jfp' q 'r']y'[f^'q''r''\y''\[fpqr]y) 

where the quantities K are isoscalar factors and S are CG coefficients. In particular the last 
factor is the CG of S3. We use the same notations as in Ref. |L7}] : [f p ] corresponds to the 
partition of S5 obtained after removal of the particle 6, [f pq ] to the partition of S4 obtained 
after removal of the particle 5, etc. 
The 3-body cfp is defined as: 

mfyq , rV , Vq''r''\[f}pqr) = K {[f]p'[f"]p"\[f]p) K {[f'Mf^VM) ™ 

x^([/; v ]r'[^ v ^"i[^]0 

The values of the K3 can then be calculated by using the corresponding tables of Ref. 



|T7| . They are listed in Tables |TT| and [TV]. They give respectively the cfp relevant for the 
decomposition of the [222] oc and [21 4 ] OC ; state. In the calculation of the expectation values, 
the CG of S3 are not necessary, as they are added up in the orthogonality relation. 
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TABLES 

TABLE I. Examples of the matrix elements (ip n \d abc F? F b Ff: \tp n ) for a few values of the indices 
(ijk) at fixed (abc) = (146). The successive columns correspond to the states ipi,...,ip5; the last 
column gives the sum over the five states. 



(ijk) 


ipi 


^2 


V>3 




V>5 


5 

E 

n=l 


(123) 


1/48 


-1/384 


-1/384 


-1/384 


-1/384 


1/96 


(145) 





-1/768 


1/768 


1/768 


7/768 


1/96 


(124) 





7/384 


-1/384 


-1/384 


-1/384 


1/96 



TABLE II. Values of the non-vanishing matrix elements (ip2\d a c F^F 2 F^\ip2) ■ The first column 
gives the colour indices (abc) ; the second the corresponding constant d abc and the third column the 
value of the matrix element. 



(abc) 


d abc 




118 


V3/3 


-1/288 


146 


1/2 


-1/384 


157 


1/2 


-1/384 


228 


V3/3 


-1/288 


247 


-1/2 


-1/384 


256 


1/2 


-1/384 


338 


v/3/3 


-1/288 


344 


1/2 


-1/384 


355 


1/2 


-1/384 


366 


-1/2 


-1/384 


377 


-1/2 


-1/384 


448 


-V3/6 


-1/1152 


558 


-V3/6 


-1/1152 


668 


-V3/6 


-1/1152 


778 


-V3/6 


-1/1152 


888 


-V3/3 


-1/288 
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TABLE III. The 3-body cfp K 3 ([42]p'q'r'[222]p"q"r"\ [222]pqr). The rows correspond to p'q'r' 
and the columns to p" q"r" . The value of pqr is given in the upper-left corner of the table. 



pqr = 332 


332 


323 


321 


221 


V 2 / 12 






212 




^10/108 




211 




-^10/108 


^10/108 


122 




-^5/108 




121 




^5/108 


-^5/108 


112 


V 5 / 12 







pqr = 323 


332 


323 


321 


221 




^4/54 




212 


^10/108 


-^5/162 




211 


-^10/108 


-^20/162 


-^5/162 


122 


-^5/108 


-^10/162 




121 


^5/108 


-^40/162 


-^10/162 


112 




-y/ 5/108 


^5/108 




pqr = 321 


332 


323 


321 


221 






^4/54 


212 






^20/162 


211 


-^20/108 


^/l0/162 




122 






^40/162 


121 


v/5/54 


^20/162 




112 




75/54 
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TABLE IV. Same as Table [D but for K 3 ([42]pq'r'[222]p"q"r"\ [2l A ]pqr). 


pqr = 543 


332 


323 


321 


221 




/ , 
V 5 /27 




ill 


-"V 4 / 27 


/1 /ei 
-y 1/81 




211 


\A/27 


-^4/81 


Vvsi 


122 


^2/27 






121 


-^2/27 


-^8/81 


V 2 /81 


112 




^2/27 





pqr = 541 


332 


323 


321 


221 






^5/27 


212 






^4/81 


211 


\A/27 


\/2/81 




122 






^8/81 


121 


V 4 / 27 


^4/81 




112 




^4/27 





pqr = 154 


332 


323 


321 


122 






y/l/E 


121 




-<\/2/5 




112 


y/2/5 







pqr = 514 


332 


323 


321 


212 






V 5 /27 


211 




sj 10/27 




122 






-^8/135 


121 




^16/135 




112 


^4/15 
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